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Slab ASE 01 
Inputs: 

10 
0.35 
0.25 
1.82 
1.6 
0.08 
100000 


Outputs: 
0.08 
-21.9501 


slab length (cm) 
slab height (cm) 
slab thickness (cm) 
slab refractive index 
parasitic coating index 
specific gain (nepers/cm) 
number of rays to launch 



maximum gain (nepers/cm) 
minimum gain (nepers/cm) 
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Dim GainDistribution(lOOO) 
Const' pi As Double = 3.141592654 

Sub Main() 

' Main Macro 

1 Macro recorded >y Raymond J. Beach 

• Keyboard Shortcut: Ctrl+u 


•Get input parameters 

Worksheets ("sheetl") .Select 

Range ("length") -Select: SlabLength - ActiveCell .Value 
Range ("height") .Select: SlabHeight = ActiveCell .Value 
Range ("thickness") .Select: SlabThickness = ActiveCell. Value 
Range ("slabindex") .Select: Slablndex = ActiveCell .Value 
Range ("coatingindex") .Select: Coatinglndex = ActiveCell. Value 
Range ( "specif icgain") .Select: SpecificGain = ActiveCell .Value 
Range ("numberof rays") .Select: NumberOfRays = ActiveCell. Value 

■Define other parameters 

NumberOfParasiticDirections = 0 
Nbins = 100 

MaxGain = SpecificGain 

Range ("maximumgain K ) .Select: ActiveCell. Value = MaxGain 
Relativelndex = Slablndex / Coatinglndex 
If SlabHeight < SlabThickness Then 

MinGain = 2 * Log ( (Relativelndex - 1) / (Relativelndex + 1) ) / SlabHeight 

Else 

MinGain = 2 * Log ( (Relativelndex - 1) / (Relativelndex + 1) ) / SlabThickness 
End If 

Range ("minimumgain") .Select: ActiveCell. Value = MinGain 

'Initialize the random number generator 
Randomize 

•Start the launch cycle 
For i = 1 To NumberOfRays 

'Define a random launch direction in {+,+,+) quadrant using direction cosines to define the dire 
tion 

Phi = (pi / 2) * Rnd 

Theta = (pi / 2) * Rnd 
*x is the slab height direction 
•y is the slab thickness direction 
1 z is the slab length direction 

cx = Sin(Theta) * Cos (Phi) 'direction cos in x-direction 

cy = Sin (Theta) * Sin (Phi) 'direction cos in y-direction 

cz = Cos (Theta) 'direction cos in z-direction 

•Define unpolarized Fresnel reflection coefficients for three different planes that generate im 
ge space 

* x-plane calculation 

Theta 1 = ArcCos(cx) 

Temp = Slablndex * Sin(Thetal) / Coatinglndex 
If Abs(Temp) > 1 Then 
Refx = 1 

Else 

Theta2 = ArcSin(Temp) 

Refx = ((Sin(Thetal - Theta2) / Sin(Thetal + Theta2) ) A 2 + (Tan(Thetal - Theta2) / Tan 
Thetal + Theta2) ) A 2) / 2 
End If 

• y-plane Calculation 
Thetal = ArcCos(cy) 

Temp = Slablndex * Sin (Thetal) / Coatinglndex 
If Abs(Temp) > 1 Then 


Refy = 1 

Else 

* Theta2 = ArcSin (Temp) 
Refy = ((Sin(Thetal - Theta2) / Sin(Thetal + Theta2) ) A 2 + (Tan(Thetal - Theta2) / Tar 
Thetal + Theta2)) A 2) / 2 
End If 

1 z-plane calculation 
Thetal = ArcCos(cz) 
Temp = Slablndex * Sin (Thetal) 
If Abs(Temp) > 1 Then 
Refz = 1 

Else 

Theta2 = ArcSin (Temp) 

Refz = ((Sin(Thetal - Theta2) / Sin(Thetal + Theta2) ) A 2 + (Tan(Thetal - Theta2) / Tar 
Thetal + Theta2)) A 2) / 2 
End If 

■Calculate the loss per cm in nepers/cm due to x, y, and z reflections 
Nepersx = cx * Log(Refx) / SlabHeight 
Nepersy = cy * Log (Refy) / SlabThickness 
Nepersz = cz * Log (Refz) / SlabLength 

•Calculate the net gain-loss in nepers/cm seen by this ray 
Nepers = SpecificGain + Nepersx + Nepersy + Nepersz 

'Bin this launch 

BinNumber = Nbins * (Nepers - MinGain) / (MaxGain - MinGain) 
If BinNumber < 0 Then BinNumber = 0 

GainDistribution (BinNumber) = GainDistribution (BinNumber) + 1 


If Nepers > 0 Then 
Beep 

NumberOfParasiticDirections = 
Worksheets ( t, sheet2 n ) .Cells (1 
Worksheets ( n sheet2 n ) . Cells ( 1 
Worksheets Csheet2") .Cells (1 
Worksheets ( "sheet2 n ) . Cells ( 1 
Worksheets ( n sheet2 n ) .Cells (1 
Worksheets ( n sheet2 n ) .Cells (1 
Check = Sqr(cx A 2 + cy A 2 + cz 
End If 


NumberOfParasiticDirections + 


NumberOfParasiticDirections, 
NumberOfParasiticDirections, 
NumberOfParasiticDirections, 
NumberOfParasiticDirections, 
NumberOfParasiticDirections, 
NumberOfParasiticDirections, 
2) 


1 

1) 

2) 

3) 

4) 

5) 

6) 


.Value 
.Value 
.Value 
.Value 
.Value 
.Value 


cx 
cy 
cz 

Re fx 
Refy 
Refz 


Next i 


End Sub 

^Function ArcCos(C) 
'Returns the Arc Cos of C. 

If C = 0 Then 

ArcCos = pi / 2 

Else 

ArcCos = Atn(Sqr(l - C A 2) / C) 
End If 

2nd Function 

inction ArcSin(S) 
'Returns the Arc Sin of S 

If S = 1 Then 

ArcSin = pi / 2 

Else 

ArcSin = Atn(S / Sqr(l - S A 2) ) 
End If 

End Function 


